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SUMMARY 


Von Harman' s equations for thin plates with large deflections are 
solved for the special cases of rectangular plates having ratios of length 
to width of 1.5 and 2 and loaded by uniform normal pressure. The boundary 
conditions are such a3 to approximate panels with riveted edges under 
normal pressure greater than that of the surrounding panels. Center 
deflections, membrane stresses, and extreme— fiber bending stresses are 
given as functions of the pressure for center deflections up to twice the 
•thickness of the plate. For small deflections the results are consistent 
with those given by Timoshenko. 


INTRODUCTION 


A general numerical method for solving Von Kannan s equations for 
thin plates with large deflections was developed in reference 1. A 
square plate loaded by uniform normal pressure with simply supported edges 
was studied. The boundary conditions approximate the case when a riveted 
sheet-stringer panel is under normal pressure greater than that of the 
surrounding ones. 

It was subsequently decided to extend the investigation to rectangular 
plates of various ratios of length to width. Two special cases are studied- 
in this report; namely, rectangular plates for which the ratio of length 
to wid-th has the value of 1.5 or 2. (For rectangular plates having a 
length equal to or greater than twice their width, experimental evidence 
(reference 2) indicates that they can practically be regarded as infinitely 
long.) Center deflections, membrane stresses, and extreme-fiber bending 
stresses are given as functions of the pressure for center deflections up 
to twice the thickness of the plate. For small deflections, the results 
are consistent with those obtained by Timoshenko (reference 3). 

This work was conducted at Brown University under the sponsorship 
and with the financial assistance of the National Advisory Committee for 
Aeronautics . 

The author is indebted to Professor W. Prager for his kind interest. 
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SYMBOLS 


a, b 


length and width of plate; a* shorter side of 
rectangular plate 


h 


thickness of plate 


x, y, z coordinates of a point in plate 

u, v horizontal, displacements in x— and y-directions of 

points in middle surface; nondimensional forms are 
ua/h 2 and va/h 2 , respectively 

w deflection of middle surface out of its initial plane; 

nondimensional form is w/h 


P 


E, u 
D 




by 2 


normal load on plate per unit area; nondimensional. form 
is pa^/Eh^ 

Young’s modulus and Poisson's ratio, respectively 


flexural rigidity of plate 


Eh- 


12(1 - u 2 ) 


^ 5 4 a 4 

V = &c 4 + 5x 2 Sy 2 + Sy 4 



'X > 



T " 

xy 


e r e r -y f 

e x * e y > ogr 

_ II _ »! _ t» 

x 9 € y * ?xy 


F 


membrane stresses in middle surface; nondimensional 
forms are c x 'a 2 ^Sh 2 , a^'a^Eh 2 , and T^'a^/lSh 2 , 
respectively 

extreme— fiber bending and shearing stresses; nondimensional 
forms are a^'a^/Eh^, c "a 2 yish 2 , a-nd* r^'a^/sh 2 .* 
respectively 

membrane strains in middle surface; nondimensional forms 
are € 'a 2 /h 2 , e ? a 2 /h 2 , and 7 'a 2 /h 2 , respectively 

extreme— fiber bending and shearing strains; nondimensional 
forms are e x "a 2 /h 2 , € y "a 2 /h 2 , y^'a^jh 2 , respectively 

stress function; nondimensional form is F^Sh 2 
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A, A 2 , . . . A n first-, second-, . . . and nth-order differences, 

respectively 




first-order differences in x— and y-directions, 
respectively 


FUNDAMENTAL EQUATIONS 


The deformation of a thin plate, the deflections of which are large 
in comparison with its thickness hut are still small as compared with the 
other dimensions, is governed hy Von Karman’s equations: 


Str 2 dy 2 




E 


& 2 w S 2 w ] 

s^x dyy cbc^ by 2 


+ ^ aSr + = £ + h a 2 v , a s w 

ox 1 *' dx 2 dy 2 oy^ D D \dx 2 5y 2 5y 2 Sx 2 



bx by bx by 


where D = 


Eh : 


12(1 - u 2 ) 


The medianr-fiber stresses are 


t _ 


cH? 

Sy 2 

d 2 F 

bx 2 




( 2 ) 


T f 

xy 


dx dy 


and the median-fiber strains are 
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«x' = 


€ * 

y 1 


xy 


1 fdb 

\c? 

>/£l 

®W 2 


- n 




av 

ax 2 

aV 




l _ 


2(1 + \x) b 2 ? 

E bx. by 


(3) 


The extreme— fiber bending and shearing stresses are 

„ Eh 


xy 


2(1 - 

n 2 ) 

\^Sx 2 

Eh 


/a 2 v 

2(1 - 

, 2 ) 

\5y 2 

Eh 


b 2 w 

2(1 + 

n) 

c)x ^y 


+ M- 


+ M- 


•v 2 i 

d w 


by^j 


bx c 


> 


(4) 


For a riveted panel under normal pressure greater than that of the 
surrounding panels, the boundary conditions are formulated in reference 1 
and are as follows: 


along y = 


- | and 


w = 0 


by 2 


- \i 


b 2 w 
dy 2 

c* 2 


= 0 


= 0 


r y i 

J y=o 


d S F 


d^F 

cfec 2 5y 2 


- n 



> 


(5) 


dy = 0 
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along any line x = Constant . 



along any line y = Constant. 


These expressions can "be made nondimens ional "by writing 

2 / v2 

irdr) . h and «(r) in place, of F, w, p, a, 


—p f, and in place, of F, w, p, a, x, y, and 

respectively, where a is the smaller side of the rectangular plate. 
(These latter symbols are used to effect a reduction in the amount of 
writing involved in the equations.) The resulting differential equations 
can then "be transformed into finite-difference equations. In terms of 
finite differences, the differential equations (l) are replaced hy the 
following equations: 


A V + 2A % + A *¥ = 
X xy y 


A w 
< xy , 


A/-VA fw 


A v w + 2A ' 
x xy 


2 k 

w + A w = 


-- 10.8(AZ) 1 'p + 10.8(1 l.A^v + A y 2 F ofV 


? Ay F V^. 



6 


NACA TN No. 1462 


p 

where A l - Ax = £y and \i is taken to "be 0.1, which value is character- 
istic for aluminum alloys. 

In terms of finite differences, the boundary conditions are 


w m,k = 0 


(V v ) m , k - 0 

(V* - ^ F )m,k = ° 


> 


( 8 ) 


along y = 


i and. 


k— 1 




= 0 


where n = k denotes points along the edges y - £ — , n = 0 denotes 

points along the center lines y = 0, and m = i denotes any point along 
the line x = Constant in the plate. Similarly, along x = + the 

boundary conditions are 


w k,n = 0 


(a 2 wY - 0 

\ X /k,n 

(A 2f - nA =0 

\x y /k,n 




(9) 


and 



where m = k denotes points along the edges x = £ m = 0 denotes 

points along the center line x = 0 and n = i denotes any point along 
the line y = Constant in the plate. 
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The normal median— fiber stresses are 

V* 


t 


! 


(AZ) ; 


y 


2 r 

v f r 

(AZ) 2 


and the extreme— fiber bending stresses are 


(10) 


2(1 - n a )(Al) 2 


2(1 - h 2 )(AZ) 2 


( A x 2v + 

(a^^w + mA^w) 


V 


( 11 ) 


RECTANGULAR PLATE WITH LENGTH-WIDTH RATIO OF 1.5 
UNDER UNIFORM NORMAL PRESSURE 


The plate is divided into 24 square meshes. (See fig. 1.) The 
points 2 f , 5 f , 8* , 9 f , 10* , 11 1 , and 11" are fictitious points outside 
the plate in order to make possible a better approximation to the boundary- 
conditions. Since the plate is symmetrical with respect to the center 
lines, only one— quarter of the plate needs to be considered. In order to 
get satisfactory results in the subsequent computations, it is convenient 
to use a number of figures beyond those normally considered justifiable 
in view of the precision of the basic data. 
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With p 2 = 0.1 or p = 0. 316228 9 the compatibility equations 
lecome 

20F q - l6F 1 + 2F 2 - l6F 3 4- 8F^ + 2Fg - K Q 

-8F q + 21F 1 - 8F 2 + 4F 3 - l6F^ + 4F^ + 2F^ + Fg, = K ± 

-8F q + 4F-L + 2 IF ^ - l6F^ + 2F^ - 8F5 + 4 F t 4- F^ = 


2F 


0 


SF-l 4 - 2F 2 - 8F 3 + 22F^ ~ SF^ + 2Fg - 8 F t + 2Fq + F 1Q 4 - F^, = 


y 


Fq — 8F 3 + 4F ^ + 20F^ — l6Fy + 2Fg — 8F^ + 4 F^q 4 - F^, — E^ 

F 1 + 2F 3 - 8F 4 + 2F 5 - 8F 6 4- 21Fy - 8F 8 + 2F 9 - 8F 10 + 2 F i;l 4- Fq 8 + F 1Q 


t = K 


7 

( 12 ) 




2 ' 

where K Q , K^, Ky K and are equal to (A^w^ 2 - A^ wAy £ 

at points 0, 1, 3 * 4, 6, and ly respectively. 

The equilibrium equations are 

20 wq — l6w-^ — l6w 3 + 2 w 2 4- 8w^ + 2wg = p* + 10 .8 ja* 0(^3 — 2w o) 

+ P'o(^i - 2w 0 ) " 2r 'o6 U + w 0 - W 1 - W 3 )] 

21w 1 - 8w q - 8w 2 - l6w^ + 4-w, + 4-w^ + 2wy + w 2 , = p' + 10.8 jct' 1 (2w^ - 2w-^ 

+ p'Vwg - 2 w x + w 0 ) - 27' 1 (w 5 + w x - w 2 - v*)] 

21w 3 — 8 wq — l6yj^ — 8wg + 4v^ + + 4 v^. + = P* + 10 .8 [cl* “ 2w^ + ™0j. 

+ P' 3 (2w^ - 2wA - 27' 3 (w 7 + w 3 - w 4 - wg)J 


22w, + 2w Q - 8-w x + 2w 2 - 8 w 3 - 8w ? - 8w 7 + 2wg + 2wg + w 10 + w^' = p' 


Y 


+ 


10.8 


a 1 i±(wrj — 2w^ + w 1 


_) + p \( w 5 - ^4 + w 3) ~ 2r \6 


Vn + V), Wcr W 


8 


4 ~ 5 


'?)] 


20 wg + Wq — &w 3 + 4 vp — 16 wy — 8 -w^ + 2 wq + 4 w 10 + w^, - p’ 

+ 10.8ja* 6 (w 3 - 2 w 6 + + P'g(2wY “ ~ 27 ’6(^io + " w 7 ” w 9)] 

21v^ 4- w^ + 2 w 3 — 8w^ 4- 2 w 9 — 8w^ — 8 wq + 2w 9 — 8-w^q 4- 2w-^^ + W Q» + w x0 T = ^ 


10.8 


a 


’76 10 _ 2w T + w 0 + p, 76 8 " ^7 + w 6) “ 27, 7( Wl1 + w 7 - w 8 _ w lo) 


( 13 ) 
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whore p' = 12(1 - n 2 )(Al)S> = 0.042l875p since Al = and a', p«, 

and 7 * are Ay?F, and A^J 1 at the points Indicated fcy the 

subscripts, respectively. 

The conditions for zero edge displacements are ^ 

( 2 F 3 - 2F q ) - - 2F 0 ) + (bT k - - n(a ? 2 - 4F X + 2F 0 ) + (2F ? -2F^ 

- ^( F l " 2F 2 + F 2 t) - S 1 

(F 0 - 2F 3 + F 6 ) - |i (2F U - 2F 3 ) + (2F ? - 4F U + 2FJ - ^(SF^ - 4F U + 


( F 2 - 2F 5 + Fg) - n(*\ - 2F ? + F 5 ,^ - _ 2 


s. 


C f 3 - ^6 + F 9 ) - ^(^7 - &&) + (»4 -^ 7 + ^ lo ) - *(*6 - 4f 7 + ^ e ) 

+ (f 5 - 2F 8 + r u ) - ^ 8 , - 2 Fq + f 7 ) - s 3 

(fl - »o) “^3 - ®o) + (^4 - 4F 3 ) " ^ (fo " 1|F 3 + »6) + ( 4F 7 - 4F e) 

1 +F^ + 2F 

(f 2 - 2F 1 + F Q ) - n(2F^ - 2F^ + (2F ? - 4F^ + 2F^ - - 4F^ + 2F.J) 

+ (2Fg - 4F ? + 2Fg)- h( 2F^ - 4F 7 + 2F 10 ) + (F n - 2F 10 + F g ) 

2 Ft + F«, ^ - S 


> 


- h( 2 F 3 - kr 6 + 2F 9 ) + (2F 10 - SFg) - h(f 6 - 


2F rt + F 


9*) = S 4 


H ^ 7 2 F 1q + ? 10 »^ - 


(14) 


J 


where 


S, = 


= 


a. - 


Sc — 


( W 2 " W l) + ( W 1 " W o) : 

- w 4 ) 2 + ( v 4 - w 3 ) 

(^8 “ w 7) 2 + ( w 7 “ w 6 ) ! 

^9 " W S) S + ( w 6 " v 3) 2 + ( w 3 " w o) 2 
( v 10 ~ w t) 2 + ( v 7 - w 4) 2 + (wi, - w-^ 2 
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The “boundary conditions are 


v 2 = 0, w 5 = 0, Wg = 0, v g = 0, v 10 = 0, w n = 0 

V “ ^2 + V 1 = ° 

w^ 1 - 2w^ + Vj^ = 0 

wg f - 2 Wq + = 0 

V “ ^9 + w 6 = 0 
w 10- " *io + w 7 = ° 

F 2 , - 2Fg + - H ^2F 5 - 2F 2 ^ - 0 

V - 2F 5 + F^ - ^(Fg - 2F 5 + Fg) = 0 

Fg, - SFq + F ? - u(F u - 2Fg + F 5 ) = 0 

F 9' - ^9 + F 6 - ^(flO ~ *9) = ° 

F lo» " *10* + F 7 ~ K F U " ®10 + F 9) = ° 

Now the boundary— value problem determines the values of w uniquely and 
the values of F to within an unknown constant. Since the actual value 
of the constant is irrelevant, it may be defined by letting F^ = 0. 

On combining the boundary conditions with equations (12), (13), and (l4) 
the final equations are 
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20F Q - 16f x + 2F 2 - l6F 3 + 8F^ + 2F 6 = Kq 

- 8F 0 + 20F 1 - 6.632456f 2 + 4f 3 - i 6 f^ + 4.632456f^ + 2f 7 = K ± 

- 8 f 0 + 4f x + 21 f 3 - 16 f^ + 2 F 5 - 8 f 6 + 4f t + f 9 = k 3 

2F 0 - 8 f i + 2.316228 f 2 - 8 f 3 + 2 U^ - 6.632456f 5 + 2F£ - 8 f 7 + 2.316228FQ 

+ *10 - K 4 

F 0 - 8 f 3 + 4f^ + 19Fg - ^F^ + 2 Fq - 6.632456F g + k.6l2h%J 1Q = Kg 

Fi + 2F, - 8 f^ + 2.3162 2 8F 5 - SFg + lSF^ - 6.632456F3 + 2.3l6228F g 

- 6.632456f 10 = K ? 

- 2Fq - 3.367^F x - 2.432456f 2 + 2F, + 4f^ + 1.8f 5 = S-,^ 

F 0 + 2F X + 0.9F 2 - 2F, - 3.3675 ^^ - 2A32^56f 5 + Fg + 2F ? + 0.9 Fq = Sg 

F 3 + 2F^ + 0.9F 5 - 2Fg - 3.36754^ - 2.432lf56Fg -U F 9 + 2F 10 + F lx = S 

- 2F 0 + aFj^ - 4 f 3 + 4f^ - 3 . 3675^Fg + 4F ? - 2.432i*56F g + 1.8 f 10 = s g 

F 0 - 2F X +.Fg + 2F 3 - 4f^ + 2F 5 + 2Fg - 3 . 3675^F 7 + 2Fq + 0.9F g 

- 2.432^56f 10 - s 5 


and 
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{20 + 21.6(0* 0 + 3* 0 + r'q)] w 0 - [16 + 2i.6(p* 0 + r' 0 )j w i 

-I16 + 21.6 (a' 0 + 7' 0 jj w 3 + (8 + 21.67'^w^ + 2w g = p’ 

- (8 + 10.83’^ v 0 + [20 + 21.6(0^ + 3*1 + r*i)J v x + W 3 


- jl6 + 21.6( a . 0 + 7’gjj w u + 2W ? = P* 

- (8 + 10.8a' 3 )w 0 + + [21 + 21.6(i’ 3 + 3' 3 + 7* 3 )J "3 


[16 + 21.6(3’ 3 + 7 * 3 )] v u - [ 

8 + 10.8^a f 3 + 27 t 3 Jj 

(b + 21*6y * p* 

— - 


2w q — (8 + 10.8a'^ Wl - (& + 10 . 83 '^w 3 + j21 + 21.6(a'^ + 3'^ 
+ 2vg - [j3 + 10.8(a' 4 + 27 * = p' 

„0 _ (8 + 10.8a' 6 )w 3 + + [19 + 21.6(a' 6 + 3» 6 + 7’gjj « 6 

- [l6 + 21.6^3’g + 7*g)l * 7 = p' 

W 1 + ^3 ~ (® + lO.Sa’^W^ ~ ^8 + 10.83* ■yt V ^ 







By following the procedures outlined in reference 1, equations (15) and (l6) 
can be solved simultaneously for w and F by the method of successive 
appr oximat i ons . 


RECTANGULAR PLATS WITH LENGTH-WIDTH RATIO OF 2 
UNDER UNIFORM NORMAL PRESSURE 


The plate is first divided into 8 square meshes and then into 32 square 
meshes. (See figs. 2 and 3, respectively.) On referring first to figure 2, 
points l f , 3*, 4 f , and 5 " are fictitious points outside the plate in 
order to give a better approximation to the boundary conditions. Consider 
one— quarter of the plate; the compatibility equations are 

20F Q - 16F X - 16F 2 + 8F 3 + 2F^ 4- 2F 1 * * K Q 

20F 2 - 8E q - 16 F 3 - 8F 4 + 4F X + AF^ + F 2 + 2F 3 + =1^ 
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|M 2 - v^v 2 ! 


at points 0 and 1, 


vhere aid are equal to 

respectively. 

The equilibrium equations are 

20w q — l6w 1 — l6w^ + 8w^ + 2w^ 1 + 2w^ » p’ +10.8 ja f 0 ^2w 2 — 2Wq^ 

+ P'o^l - ^o) - 2 ?’o( w 3 + w 0 - w l - w 2 )j 

20w 2 — 8wq - l6w^ - 8w^ + 4v 1 + 4w^ + w 2 + 2w^ f + v^r 

“ P* + 10 - 8 p 2 (-o “ ^ + v k) + 01 2 “ ^2) 

- 2 7 ' 2 (v 5 + w 2 - w 3 - wj| 

where p 1 = 12(1 — u 2 )( Al)^ = O. 675 P since Al ~ and a*, 0 f , and 

y i are A^/^F, A^ 2 F, and ^ F the points indicated hy the subscripts, 

respectively. 


> (18) 


The conditions for zero edge displacements are 
(2F 2 - 3?o) - n^F x - 2F Q ) + (2F 3 - 2F^ - h(F q - = v Q 2 

(F 0 - 2F 2 + F^ - 11 (2F 3 - 2F 2 ) + ( f x - 2F 3 + jJJ - h(F 2 - 2F 3 + F^ 
“ ^0) " ~ ^0) " ^( F 0 " 2F 2 + F l) + 2 ^ F 3 - ^2) 

- < F 2 - 2\ + \.) + (® 5 - ^4) " ( W 2 - W o) 2 + w 2 2 


n 


= M. 


>(19 
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The boundary conditions are 


"w-i = 0, Wo - 0, w ( = 0, = 0 


w 1 
W 1 

+ w Q = 

= 0 

w . 

+ w^ : 

= 0 

3 

2 


w. 

+ w 

= 0 

4 * 

2 


F n . 

- 2F_ 

+ F 

1* 

1 

( 

F. , 

- 2 F 0 

+ F, 

3 

3 


F v 

- & k 

+ F. 


0 


The problem can now be solved uniquely for the values of w and the values 
of F to within an unknown constant. As the actual value of the constant 
is irrelevant, it may be defined by letting F = 0. 

Combined with the boundary conditions, equations (17), (l8), and (19) 
are then 


18F q - 13.264912*1 - loF p + 9.264912F 3 + 2F^ = K Q 

-8F 0 + 4 . 632456 F, + l8F 2 - 13.264912F 3 - 6.632456 f^ = K ± 

-1.367544F 0 - 2.532456F x + 2F p + 1.8F = w Q 2 

F Q + 0.9F-L - 1.367544F - 2.532456F 3 + F^ w p 2 

- 2F 0 + 2F X - 3 • 367544F2 + 4F 3 - 2.532456F^ = (w 2 - w Q ^ + w 2 2 


> 


(20) 


jx8 + 21 .6{a' 0 + p* 0 + 7'o)] w O - [16 + 21.6(a’ 0 + 7*q)] w 2 p * \ 

-(8 + 10 . 8 a * 2 ^ w 0 + fl 8 + 21.6 ( a * g + p * 2 + 7’ 2 )] w 2 = P ' 


( 21 ) 
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Equations (20) and (21) can now be solved simultaneously for w and F 
by the method of successive approximations. 

With the values of w thus computed as a first approximation, the 
case wherein the plate is divided into 32 square meshe3 is, now to he 
studied. On referring to figure 3, points 2', 5*, 8', 11’, 12’, 13' 

14* , and 14" are again fictitious points outside the plate in order 
to give a letter approximation to the boundary conditions. Consider 
one— quarter of the plate; the compatibility equations are 


20E Q - 16F 1 + 2F 2 - 16F 3 + 8F^ + 2Fg 
-8 F q + 21F 1 - 8F 2 + 4F 3 - l6F 4 +'4F 
-8F 0 + 4F.-J + 21F 3 - l6F^ + 2F 5 - 8F 6 
2 F q - SF X + 2F 2 - 8F 3 + 22F ]f - BF^ + 

F Q - 8F 3 + 4F^ + 20F 6 - 16SF + 2F 8 - 
h + 2F 3 - .+ 2F_ - 8Fg + 21F 7 - 8F 8 + 2F 9 - + 2F U + F^ + Fg, = K ? 

F 3 _ 8F 6 + UF ? + 20F 9 - 16F 10 + 2P U - 8F 12 + UF 13 + F^. = K g 

F * + 2F 6 - a 7 + ^8 ~ *9 + 21F 10 - *11 + ^IS ■- *13 + ' 2F n + hi* + V = ho 

( 22 ) 

where K = Kvl 2 - s\ at the points indicated by the subscripts. 


Kq 


+ 2F + F 0 t — K 

( 1 


+ 4F„, + F^ = Ko 
( 9 3 


2Fg - 8F 7 + 2Fg + F 1q + r = K u 


*9 + UF 10 + F 12 ’ K 6 


t 

/ 
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The equilibrium equations are 

20w Q — l6wj^ — l6w^ + 8w^ + 2 w 2 + 2w^ = p f + 10.8 — ^*o) 

+ P' 0 (ar 1 - 2* 0 ) - 2r' 0 (v^ - »! - w 3 + w Q ) 

— 8w q + 21w^ — 8wg + Uw^ — l6w^ + 4-w, + 2w^ + w^, = p' + 10 . 8 la' ^ — 2w^ 

+ P 'l( v 2 - ^1 + w 0 ) - 2 ’ , 'l( w 5 - w 2 ~ w 4 + w x ) 

— 8w^ + Av-J_ + 21w^ — l6w^ + 2w,_ — 8w^ + ~ p f + 10.8|a f ^^r 0 — 2w^ + 



/ 


(K - 

>- 27 ’ 3 (^ 7 - 

- W 4 " w 6 + w 3 ) 


2^ 0 — 8w 1 + 2w^ — 8w^ + 22w^ — Qw^ + 2w^ — 8w^ + 2Wq + 

+ v 5 , = p' + 10.8ja«^^r 7 - 2w u + w^J + 3' u (w 5 - 2w^ + 
- 2 A( w 8 - w 5 - w ? + w 4 jj 


W 


0 


8w^ + + 20w^ — l6w^ + 2 w q — 8w Q + W in + w 




8 


10 12 


= p' + 10.8^’ 6 (w 3 - 2w 6 + + P' 6 (2v ? - SVg) 

10 - v 7 " w 9 + w 6)] 

w + 2w — 8w. + 2w_ — 8w . + 2Tw - 8wn + 2w^ — 8w 

1 3 4 5 6 789 10 

• 8 f 7 (% - ^7 + W lo) 

+ p ' 7 ^8 - ^7 + w e) - 27 ’ 7 ( W 11 - w 8 - W 10 + w 7 )| 

8wv + 4w + 20w — l6w_, _ + 2w „ — 8w + W, ^ + v 


+ 2w + Wt _ + w 0f = p f + 10 

11 13 8 f 


w 


3 '"'6 

+ 10 


7 


10 


11 


12 


13 12 f 


•T' 9 ( W 6 “ ^9 + W 12) - ^' 9^13 - w io - W 12 + w 9 j| 

\ + *6 - ^7 + ^8 - % + 21w io - ^11 + ^12 - ^13 + + w ii< + w 


13' 


= p' + 10.8 


a'i 0 (v 7 - 2 Wl0 + Wl? ) + P’iQ^u - 2w i0 + w 9 ) 

“ 2r ' 10 (^ 1 ^ ” W 11 “ W 13 + w lo) 


J 


( 23 ) 
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where p* - 12^1 — M.^)(A2)^p - O.Ol*2l875p since A l - jj^, and a* , 3*, 

and y 1 are A^% f A %*, and A^^F at the points indicated by the 

, y 

subscripts, respectively. 

The conditions for zero edge displacements are 


2F 3 - 2F 0 - n(2F x - 2F^ + ItF^ - ^ - 2 h(F 2 - S^ 1 + + 2F ? - 2F 2 

- *(* 2 , ~ » 2 + F 1 ) = S 1 

F 0 " * 3 + F 6 " ~ *3) + *1 " 4F 4 + *7 " K F 3 " ^ 1+ + F 5 ) 

+ (F 2 - 2F 5 + F 0 ^ - h(F 5 , - 2F 5 + F^) = S 2 

F 3 _ 2F 6 + F 9 - u( 2F 7 - SFg) + 2F t - 4F ? + 2F 1Q - 2u(Fg - 2F ? + F g ) 

+ ( f 5 " ^8 + F ll) " K F 7 ' ?J 8 + P 8«) = S 3 

F 6 - 2F 9 + F 12 - ^ (*10 - *9) + 251 7 - *10 + *13 - 2 ^( F 9 - *10 + F ll) 

+ F 8 - *11 + F 14 - * ( F 10 - *11 + F u.) = S u 
^1 - *0 - ^ (*3 - *0) + 4 \ - '*3 - K F 0 - *3 + F e) + 4F 7 “ 4F 6 
- ^( f 3 - *6 + F 9 ) + 4F io " 4F 9 " 2tl ( F 6 “ *9 + F is) + *13 " *12 

~ ^9 “ *12 + F 12») = S 5 

F 2 — 2F-^ + F 0 — 4x( 2F 4 - 2F^ + 2F ? - 4F^ + 2F 3 - 2^ - 2F^ + F ? ) + 2Fg 
+ 2F 6 - ax (f 4 - 2F 7 + F 10 ) + 2F U - + 2F ? - ^ - 2F 10 + F^) 

+ F 14 " *13 + F 12 ~ ^( F 10 - *13 + F 13') “ S 6 

vhere \ = V (vi.n “ w m,n) S for 1 = lj 2> 3 5 4 

m=0 


and 



^in,n+l “ w m,n) 


for J = 5, 6 
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The boundary conditions are 

w 2 = 0, v 5 = 0, w 8 = 0, v u = 0, v 12 = 0, w 13 = 0, v l4 = 0 
w 2 , - 2w s + w x = 0 

Vj-, — 2w + v = 0 

5 5 4 

v 8> - ^8 + w 7 = ° 


w ll» 


+ 

II 

O 

: 0 



H*= 

ro 

~ ^XS 

+ 

V = 

9 

0 



W 13 ’ 

-^3 

+ 

w xo = 

: 0 



F l- 

2F 2 + 

F 2« 

-•M 

^2F 

- ?f 2 ) 

= 0 

F ^- 

2F 5 + 

F 5 ’ 

-•*( 


• 2F 5 + 

0 

II 

F 7 - 

2F 8 + 

F 8* 

-n( 


•2F 8 + 

F x^) -° 


F 10 " + F xx» " •* ( F 8 ~ + F ik) = 0 

F 9 " ^xa + F i2' - • , ( 2p 13 " ^ 12 ) = 0 
F xo - ®13 + f 13' " u (> - ^13 + F ia) ■ 0 

The problem can now be solved uniquely again for the values for w and the 
values of f to within an unknown constant. As the actual value of the 
constant is irrelevant, it may be defined by letting = 0. 

On combining with the boundary conditions, equations (22), ( 23 ) , and (24) 
are then 
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20F q - 16f x + 2F 2 - 16f 3 + aF^ + 2Fg = j 

- 8 f 0 + 20^ - 6.632456F2 + 4f 3 - i 6 f 4 + 4.632456f 5 + 2 f ? = k x 

- 8 f 0 + 4F-L + 21F 3 - l6F^ + 2Ftj - 8Fg + 4 f ? + Fg = K g 

^0 - ^1 + 2.316228 f 2 - 8f 3 + 21F^ - 6.632456F g + 2Fg - 8 f 7 + 2.316228F 

+ F 10 = k 4 

F 0 - 8F 3 + 4F 4 + 20Fg - 16F ? + 2F g - 8F 9 + 4F 1q + F lg = Kg 

Fi + 2F 3 - 8F^ + 2.3l6228F 5 - 8Fg + 20®^ - 6.632456Fg + 2F g - 8 F 1q 

+ 2 . 3l6228F 1 - + F., = K 
11 13 7 

F 3 “ ^6 + I+F y + 19F 9 ~ lfiF io + ~ 6.632456F i2 + 4.632456 f 13 = Kg [ 

F 4 + ^6 ~ ^7 + 2 -3i6228F 8 - 8F g + 19F 1Q - 6.6321^^ + 2.3 i6228f 12 / 

-6.632456 F = K 1q 

-2?0 - 3.367544F l - 2.432456f ? + 2F 3 + 4F^ + i.8f = s i 

F 0 + ^’l +0 -9 F 2 - 2F 3 ~ 3.36745 ^ - 2.432456f ? + Fg + 2JF + 0.9Fg = s 
F 3 + 2F 4 +G9F 5 - 2Fg - 3.36?456 f^ - 2.432456F g + F 9 + 2F 1Q + 0.9F U = S 

F 6 + 251 7 + q 9 f 8 - & 9 - 3.36?4 5 6F 10 - 2.4 3 24 5 6F 11 + v ls + 2F 13 - S u 


-2F q + aFj^ - 4F 3 + 4F^ - 4F g + 4F 7 - 3.367544F g + 4F 10 - 2.432456F 12 

+ 1 ' 8f 13 = S 5 

F 0 - 2Fi + F 2 + 2F 3 _ 4F 4 + 2F 5 + 2F g - 4F^ + 2F g + 2F g - 3-367544F 10 

+ 2 F 1X + 0 . 9 F 12 - 2.432544F = Sg 
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and 


{20 + 21 . 6 (a ' c + 3' 0 + 7 ’ 0 | w 0 


16 4- 21.6 


( e '» - ? 'o) 


w 


1 


4 - 2w> = p* 


W 1 + W 3 


— 16 4- 21.6 (a* Q 4- o) + (b + 10.87’^ 

(q + 10. 83'^ V 0 + jso + 21.6(0.^ + 3' 1 + 7*^ 

^6 + 21.6^a , 7 4 - 7*q^j v ^ + 2w = P* 

•(8 4- 10 . w q + 4 w. + jj?l 4* 21 . 6 /a* ^ 4- p * ^ + 7*^j 

~ [16 + 21.6(3*^ + 7 — 8 + 10.8(a T ^ 4- 27*3) J w 6 


4- ^ 4 21.67* ^ 4- ~ p* 

2Wq -^84 10.8a* — ^8 4- 10.8p T ^ 4 - jj51 4 - 2l.6^a*^ 4 - 3*^ 4 - 7 *^ 


4- 2w 6 


8 4- 10.8 ^L* ^ 4 - 27\J W ? 4- W 10 = p* 


w 4 


Vq “ { 8 + 10.8o.*^W- 4- 4w^ 4- |^20 4- 21.6^0.*^ 4- 3’^ 4* 7*^1 w 

16 + 21.6 (p’ 6 + r'g)] w 7 - [ 8 + 10.8^a'g + 27'gj w 9 


+ (44- 21.67*^ w 10 = p* 

w x + 2w^ -(84- 10.8a* 7 ^v 4 -(^84 10.83 * ? > 6 4- J20 4- 21.6(0,^ 4- 3* ? 4- 7* 
+ 2W9 - Js 4- 10.8^a* 7 4- 27 * 7 ) I w 10 - p* 

w 3 - (s 4 - 10.8a*^ 4- 4 -w 7 4 [19 4 21 . 6 (a * ^ + P«^ + 7*^J 
- _l6 + 2i.6(p* 0+ r'o)] w 10 « P * 

V j. + 2 w^.-r- ^8 + 10 . 8 a* 1( ^ w 7 — (s 4 10.83*^0^ 

(a* 4 3* + 


1 j 


9 


9 


+ [l9 + 21.6 (a- 1Q H 3- 10 + r - lq ) 


W 10 = p ’ 


By the method, of successive approximations, equations (29) and (26) can be 
solved simultaneously .for values of w and F. 



NACA TN No. 1^62 


21 


RESULTS AND DISCUSSION 


The "bending problem of rectangular plates under uniform normal 
pressure with simply supported edges is studied by finite-difference 
approximations. The difference equations are solved by the method of 
successive approximation. The case of a square plate is studied in 
reference 1. In the present report, the computations are extended to 
rectangular plates with ratios of length to width of 1.5 or 2. Experi- 
mental work (reference 2) indicates that a rectangular plate having a 
length equal to twice Its width may be regarded practically as an 
infinitely long plate. The maximum normal pressure calculated Is 
A 

^ a '-- — 25O, which gives the center deflections about twice the thickness 
Eh^ 

of the plate. 


The deflections at various points In the plate under different 
pressure ratios are tabulated in tables 1, 2, and 3 . The center deflec- 
tions are plotted against the normal pressure ratio in figures b and 5. 

The membrane stresses in the centers of the plates and at the centers of 
the edges are tabulated in tables b , 5, and 6 and are plotted in figures 6 
and 7. The bending and total stresses at the centers of the plates are 
tabulated in tables 7 and 8 and are plotted In figures 8 and 9- 


As indicated in reference 1, the plate should be designed for strength 
on the basis of clamped edges and for maximum deflection or washboarding 
according to the boundary conditions such as formulated in the present 
report. The finite-difference equations give the values of deflections 
with good approximations, whereas the stresses, which are obtained by a 
second-order numerical differentiation, are always less accurate. 


Brown University 

Providence, R. I. 

November b , 19^6 
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TABLE 1 

CENTER DEFLECTIONS 


L 



12.5 

25 

50 

75 

100 

150 

200 

250 


w 0 /L 


_ A _ a 

*5; az = jj- 

b _ A , ■ a 

- = 2: AZ = - 
a 2 

- = 2; Al - 
a 

0 

0 

0 

.62*4-9 

. 7322 

.6950 

.8790 

1.002 

.9^63 

1.175 

1.323 

1.241 

1.374 

1.54a 

1.440 

1.528 

1.713 

1.596 

1.769 

1.982 

1 

1.840 

1.957 

. i 

2.194 

2.033 

2.115 

2.372 

2.196 


-TIP 




















NACA TN No. 1462 


25 


TABLE 3 

DEFLECTIONS AT VARIOUS POINTS 



Ha 

BH 

W 1 

h 

^3 

h 

mm 

h 

W 7 

h 

V 

12 

h 

w 10 

h 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12.5 

.6958 

.5025 

.6707 

.4847 

.5749 

.4159 

.3585 

.2590 

25 

.9^63 

.6868 

.9201 

.6682 

.8074 

.5867 

.5184 

.3745 

50 

1.2U1 

.9055 

1.215 

.8870 

I.O87 

.7940 

.7179 

.5179 

75 

1.440 

1.054 

1.414 

1.036 

1.277 

.9354 

.8558 

, 6 l 64 

100 

1.596 

1.171 

1.569 

1.152 

1.426 

1.046 

.9646 

.6938 

150 

1.840 

1.355 

1.813 

1.335 

I.658 

1.220 

1.135 

.81 46 

200 

2.033 

1.500 

2.00k 

1.479 

1.840 

1.356 

1.269 

.9097 

250 

2.196 

1.621 

2.166 

1.600 

1.994 

1.471 

1.383 

.9889 

















TABLE 4 a 

MEMBRANE STRESSES 
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MBtsm 

a*^a 2 

yO 

■a 

rr® a 2 

a yl 

a* a 2 
x2 

mm 

Eh 2 

Eh 2 

Eh 2 

Eh 2 

Eh 2 

Eh 2 

0 

0 

0 

0 

0 

0 

1.064 

.6050 

1.072 

.3389 

.2316 

.7325 

2.109 

1.213 

2.150 

.6799 

.4775 

1.510 

3.783 

2.202 

3.903 

1.234 

.8930 

2.824 

5.176 

3.039 

5.381 

1.702 

1.251 

3.955 

6.410 

3.784 

6.694 

2.117 

1.574 

4.977 

8.605 

5.110 

9.037 

2.858 

2.154 

6.811 

10.55 

6.292 

11.13 

3.518 

2.673 

8.452 

12.34 

7.375 

13.04 

4.125 

3.153 

9.972 


Subscript 0 denotes the center of the plate; subscript 1 denotes 
the center of the sides x * i and subscript 2 denotes the 
center of the sides y = i 
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TABLE 5 a 


MEMBRANE STRESSES 


pa^ 

°'io a2 

a , _ n a 2 
yo 

a f a 2 
xl 

at yi a2 

! a *x2 a 

0* a S 
y2 a 

Eh^ 

Eh 2 

Eh 2 

! 

Eh 2 

Eh 2 

Eh 2 

Eh 2 

0 

0 

0 

0 

0 

0 

0 

12.5 

1.149 

.6330 

1.115 

.3525 

.1897 

.5998 

25 

2.132 

1.201 

2.103 

.6652 

.3751 

1.186 

50 

3.692 

2.119 

3.690 

1.167 

.6821 

2.157 

75 

4.995 

2.893 

5.024 

1.589 

.9437 

2.984 

100 

6.156 

3.585 

6.215 

1.965 

1.178 

3.726 

150 

8.218 

4 . 8l6 

8.332 

2.635 

1.597 

5.049 

200 

10.06 

5-916 

10.22 

3.233 

1.971 

6.234 

250 

11.75 

6.928 

11.96 

3.782 

2.316 

7.324 


Subscript 0 denotes the center of the plate; subscript 1 denotes 
the center of the sides x = ± and subscript 2 denotes the 
center of the sides y * i 




















TABLE 6 a 


MEMBRANE STRESSES 


KACA TN No. llt -62 


o’ ~a S 
x2 


rt* *,2 

o 

y2 


12.5 


1.293 


.5421 


2.403 1.035 

4.155 1.827 

5.612 2.496 

6.911 3.094 

9.220 4.160 

11.28 
13.16 


1.298 

2.441 


11.72 


.4106 

.7718 


4.264 1.349 


7.149 2.261 

9.570 3.026 


3.707 


.2372 


5.788 1.830 1.220 


1.534 


2.100 


2.606 


13.70 4.334 3.076 


.7500 


.4718 1.492 

.8730 2.761 


3.859 

4.852 

6.640 

8.240 

9.728 


Subscript 0 denotes the center of the plate; subscript 1 denotes 

the center of the sides x » 1 ~j and subscript 2 denotes the 

, 2 

center of the sides y m ± 










4 

pa 

a" a 
xQ 

CT "y0 a 

_.tt _ 

xO 

a, ’y0 a 

Eh 14 - 

Eh 2 

Eh 2 

Eh 2 

Eh 2 

0 

0 

0 

0 

0 

12.5 

3.413 

2.075 

4.477 

2.680 

25 

4.700 

2.761 

6.809 

3.974 

50 

6.137 

3.V76 

9.920 

5.679 

75 

7.072 

3.914 

12.25 

6.953 

100 

7.813 

4.260 

14.22 

8.044 

150 

8.909 

^.777 

17.51 

9.887 

200 

9.770 

5.177 

20.32. 

11.47 

250 

10.49 

5.519 

22.83 

12.89 


Subscript 0 denotes the center of the plate. 


















30 


NACA TN No. 1*4-62 


TABLE 8 a 

EXTREME-FIBER BENDING AND TOTAL 
STRESSES AT CENTER OF PLATE 


a 



°"xO a? 

ft” n 2 

a yO a 

°x0 

a 2 


2 

a 

pa^ 

r 

Eh* 

Eh 2 

Eh 2 

Eh 

2 ~ 

Eh 

~ 

AZ = - 
2 

cd|-=J- 

!l 

ro 

< 

AZ = § 
2 

A. O' 

AZ * 4 

AZ = - 

2 

AZ = | 

AZ = - 
2 

AZ = f - 
4 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12.5 

3.431 

3.577 

1.588 

1.532 

4.580 

4.871 

2.221 

2.074 

25 

*4-. 662 

4.761 

2.069 

1.925 

6.795 

7.164 

3.271 

2.960 

50 

6.122 

6.108 

2.625 

2.349 

9.81*4- 

10.26 

4.745 

4.176 

75 

7.116 

7.OO5 

3.003 

2.637 

12.11 

12.62 

5.896 

5.133 

100 

7.895 

7.703 

3.300 

2.86*4- 

14.05 

14.61 

6.885 

5.957 

150 

9.118 

8.792 

3.767 

3.224 

17.34 

18.01 

8.583 

7.385 

200 

10.08 

9.651 

*+.138 

3.519 

20.14 

20.93 

10.05 

8.63O 

250 

10.90 

- 

10.39 

A. A52 

3.766 

22.64 

23.55 

11.38 

9.755 

- 


a 

Subscript 0 denotes the center of the plate. 



































































Figure 8.- Extreme -fiber total and bending stresses at center of plate. 











